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It is an open challenge to analyze the crystalline color superconducting phases that may arise in 
cold dense, but not asymptotically dense, three-flavor quark matter. At present the only approxi- 
mation within which it seems possible to compare the free energies of the myriad possible crystal 
structures is the Ginzburg-Landau approximation. Here, we test this approximation on a particu- 
larly simple "crystal" structure in which there are only two condensates (us) ~ A exp(iq2 ■ r) and 
(ud) ~ A exp(iq 3 ■ r) whose position-space dependence is that of two plane waves with wave vectors 
q2 and q 3 at arbitrary angles. For this case, we are able to solve the mean-field gap equation with- 
out making a Ginzburg-Landau approximation. We find that the Ginzburg-Landau approximation 
works in the A —> limit as expected, find that it correctly predicts that A decreases with increasing 
angle between q2 and q 3 meaning that the phase with q2 || q 3 has the lowest free energy, and find 
that the Ginzburg-Landau approximation is conservative in the sense that it underestimates A at 
all values of the angle between q2 and q 3 . 

PACS numbers: 12.38.Mh,24.85.+p 



I. INTRODUCTION 

Quantum chromodynamics predicts that at densities 
that are high enough that baryons are crushed into quark 
matter, the quark matter that results features pairing be- 
tween quarks at low enough temperatures, meaning that 
it is in one of a family of possible color superconduct- 
ing phases The "laboratory" where color supercon- 
ducting quark matter is most likely to be found is the 
interior of compact stars. Except during the first few 
seconds after their birth in supernovae, these stars have 
temperatures well below the tens of MeV critical temper- 
ature expected for color superconductivity, meaning that 
if these stars feature quark matter cores, these cores will 
be color superconductors. 

The essence of color superconductivity is quark pair- 
ing, leading to Meissner effects for (some) color mag- 
netic fields, driven by the BCS mechanism which oper- 
ates whenever there are attractive interactions between 
fermions at a Fermi surface Q . The interaction between 
quarks in QCD is strong and is attractive between quarks 
that are antisymmetric in color, so we expect cold dense 
quark matter to exhibit color superconductivity. We shall 
only consider Cooper pairs whose pair wave function is 
antisymmetric in Dirac indices — the relativistic general- 
ization of zero total spin. (Other possibilities have been 
investigated [ll S H 01 l| an d found to be less favorable.) 
This in turn requires antisymmetry in flavor, meaning 
in particular that the two quarks in a Cooper pair must 
have different flavor. 



'Electronic address: massimo@lns.mit.edu 
^Electronic address: krislma@lns.mit.edu 
t Electronic address: sharma@mit.edu 



It is by now well-established that at sufficiently high 
densities, where the up, down and strange quarks can 
be treated on an equal footing and the disruptive effects 
of the strange quark mass can be neglected, quark mat- 
ter is in the color-flavor locked (CFL) phase, in which 
quarks of all three colors and all three flavors form con- 
ventional Cooper pairs with zero total momentum 0, 0] , 
and all fermionic excitations are gapped. However, even 
at the very center of a compact star the quark number 
chemical potential fj, cannot be much larger than 500 
MeV, meaning that the strange quark mass M s (which 
is density dependent, lying somewhere between its vac- 
uum current mass of about 100 MeV and constituent 
mass of about 500 MeV) cannot be neglected. Further- 
more, bulk matter, as relevant for a compact star, must 
be in weak equilibrium and must be electrically and color 
neutral H, S B S [U3| ■ All these factors work to sepa- 
rate the Fermi momenta of the three different flavors of 
quarks, and thus disfavor the cross-species BCS pairing 
that characterizes the CFL phase. If we imagine begin- 
ning at asymptotically high densities and reducing the 
density, and suppose that CFL pairing is disrupted by 
the heaviness of the strange quark before color supercon- 
ducting quark matter is superseded by baryonic matter, 
the CFL phase must be replaced by some phase of quark 
matter in which there is less, and less symmetric, pairing. 
The nature of this phase is not currently established. 

Within a spatially homogeneous ansatz, the next phase 
down in densit y is the gapless CFL (gCFL) phase |!Ull2l 
[H El E E [H [3 . In this phase, quarks of all three 
colors and all three flavors still form ordinary Cooper 
pairs, with each pair having zero total momentum, but 
there are regions of momentum space in which certain 
quarks do not succeed in pairing, and these regions are 
bounded by momenta at which certain fermionic quasi- 
particles are gapless. This variation on BCS pairing — 
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in which the same species of fermions that pair feature 
gapless quasiparticles — was first proposed for two flavor 
quark matter 0] and in an atomic physics context [2fT |. 
In all these contexts, however, the gapless paired state 
turns out in general to suffer from a "magnetic instabil- 
ity" : it can lower its energy by the formation of counter- 
propagating currents pHl22l|. In the atomic physics con- 
text, the resolution of the instability is phase separation, 
into macroscopic regions of two phases in one of which 
standard BCS pairin g oc curs and in the other of which 
no pairing occurs [231 l24l l25j . In three-flavor quark mat- 
ter, where the instability of the gCFL phase has been 
established in Refs. [^I, phase coexistence would require 
coexisting components with opposite color charges, in ad- 
dition to opposite electric charges, making it very un- 
likely that a phase separated solution can have lower 
energy than the gCFL phase 0, Furthermore, 
color superconducting phases which are less symmet- 
ric than the CFL phase but still involve only conven- 
tional BCS pairing, for example the much-studied 2SC 
phase in which only two colors of up and down quarks 
pair S_ 2(| [23] but including also many other possibil- 
ities [28 | , cannot be the resolution of the gCFL insta- 
bility [6 , [2^ . It seems likely, therefore, that a ground 
state with counter-propagating currents is required. This 
could take the form of a crystalline color superconduc- 
tor m m m m m m im e m m — the qcd 

analogue of a form of non-BCS pairing first considered 
by Larkin, Ovchinnikov, Fulde and Fcrrell |39|, Or, given 
that the CFL phase itself is likely augmented by kaon 
condensation 0, El > it could take the form of a phase 
in which a CFL kaon condensate carries a current in one 
direction balanced by a counter-propagating current in 
the opposite direction carried by gapless quark quasi- 
particles |42j |. This meson supercurrent phase has been 
shown to have a lower free energy than the gCFL phase. 
The investigation of crystalline color superconductivity 
in three-flavor QCD has only just begun [31} . Although 
such phases seem to be free from magnetic instability |38j , 
it remains to be seen whether such a phase can have a 
lower free energy than the meson current phase, mak- 
ing it a possible resolution to the gCFL instability. The 
simplest "crystal" structures do not suffice, but experi- 
ence in the two-flavor context [3(| suggests that realistic 
crystal structures constructed from more plane waves will 
prove to be qualitatively more robust. Once we know (or 
have a well-motivated conjecture for) the favored crys- 
tal structure, the challenge will then be to calculate its 
shear modulus and the pinning force it exerts on rota- 
tional vortices trying to move through it, as these are 
the microphysical inputs needed to determine whether 
observations of pulsar glitches can be used to rule out 
(or in) the presence quark matter in the crystalline color 
superconducting phase within compact stars |29j| . 

Determining the favored crystal structure(s) in the 
crystalline color superconducting phase (s) of three- flavor 
QCD requires determining the gaps and comparing the 
free energies for very many candidate structures, as there 



are even more possibilities than the many that were inves- 
tigated in the two-flavor context [3(|. As there, it seems 
likely that progress will require making a Ginzburg- 
Landau approximation. This approximation is controlled 
if A <C Ao, where A is the gap parameter of the crys- 
talline color superconducting phase itself and Ao is the 
gap parameter in the CFL phase that would occur if M s 
were zero. This approximation will be dubious precisely 
in the case of interest: if a crystalline phase exists with 
a free energy lower than that of the gCFL phase, such 
a phase will be characterized by robust pairing meaning 
that A will not be much smaller than Ao . Our purpose in 
this paper is to analyze a particularly simple one param- 
eter family of "crystal" structures in three-flavor quark 
matter, simple enough that we can do the analysis both 
with and without the Ginzburg-Landau approximation. 

We shall work throughout in a Nambu-Jona-Lasinio 
(NJL) model in which the QCD interaction between 
quarks is replaced by a point-like four-quark interaction, 
with the quantum numbers of singlc-gluon exchange, an- 
alyzed in mean field theory. This is not a controlled ap- 
proximation. However, it suffices for our purposes: be- 
cause this model has attraction in the same channels as in 
QCD, its high density phase is the CFL phase; and, the 
Fermi surface splitting effects whose qualitative conse- 
quences we wish to study can all be built into the model. 
Note that we shall assume throughout that A <C ft- 
This weak coupling assumption means that the pairing 
is dominated by modes near the Fermi surfaces, justify- 
ing the analysis of the model using High Density Effective 
Theory (HDET) OE^. Thus, we shall be comparing 
two calculations: one making the Ginzburg-Landau as- 
sumption that A <C Ao -C /i and the other allowing for 
A ~ Ao <C /i. The Ginzburg-Landau calculation can be 
extended to more complicated crystal structures than the 
simple ones that we shall analyze. We shall find that the 
Ginzburg-Landau approximation works when it should 
and that, at least for the simple crystal structure we an- 
alyze, it is conservative in the sense that when it breaks 
down it always underestimates the gap A and the conden- 
sation energy. Furthermore, we find that the Ginzburg- 
Landau approximation correctly determines which crys- 
tal structure among the one parameter family that we an- 
alyze has the largest gap and lowest free energy. We shall 
see, however, that the range of validity of the Ginzburg- 
Landau approximation does depend on the crystal struc- 
ture. 

Our paper is organized as follows. In Section[H]wc shall 
describe the one parameter family of three-flavor crys- 
talline color superconductors that we analyze. In Sec- 
tion II I II we present the results of our Ginzburg-Landau 
analysis. In Section ITVI we recast the High Density Ef- 
fective Theory slightly, as needed for our purposes, and 
redo our analysis without the Ginzburg-Landau approx- 
imation. In Section we show the numerical results of 
our analyses, make comparisons, and draw conclusions. 
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II. MODEL AND ANSATZ 



weak equilibrium are given in the absence of pairing by 



We shall analyze quark matter containing massless 
u and d quarks and s quarks with an effective mass 
M s - (Although the strange quark mass can be deter- 
mined self-consistently by solving for an (ss) conden- 
sate 0, 0, 0|, we shall leave this to future work and 
treat M s as a parameter.) The Lagrangian density de- 
scribing this system in the absence of interactions is given 
by 
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where i,j = 1,2,3 are flavor indices and a, = 1,2,3 
are color indices, where we have suppressed the Dirac in- 
dices, where = S al3 diag(0, 0, M s )ij is the mass ma- 
trix, where = dd a/3 Sij and where the quark chemical 
potential matrix is given by 
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with Q — diag(2/3, — 1/3, — l/3)jj the quark electric- 
charge matrix and T3 and Ts the Gcll-Mann matrices 
in color space. We shall quote results at quark number 
chemical potential /j, = 500 MeV throughout. 

In QCD, /i e , /i3 and [is ar e the zeroth components 
of electromagnetic and color gauge fields, and the gauge 
field dynamics ensure that they take on values such that 
the matter is neutral @, E3i satisfying 



= 



(3) 



with fl the free energy density of the system. In the NJL 
model that we shall employ, in which quarks interact via 
four-fermion interactions and there are no gauge fields, 
we introduce fi e , ^,3 and by hand, and choose them 
to satisfy the neutrality constraints Q. The assumption 
of weak equilibrium is built into the calculation via the 
fact that the only flavor-dependent chemical potential is 
/i e , ensuring for example that the chemical potentials of 
d and s quarks with the same color must be equal. Be- 
cause the strange quarks have greater mass, the equality 
of their chemical potentials implies that the s quarks have 
smaller Fermi momenta than the d quarks in the absence 
of BCS pairing. In the absence of pairing, then, because 
weak equilibrium drives the massive strange quarks to be 
less numerous than the down quarks, electrical neutrality 
requires a fi e > 0, which makes the up quarks less numer- 
ous than the down quarks and introduces some electrons 
into the system. In the absence of pairing, color neutral- 
ity is obtained with [13 = /ig = 0. 

The Fermi momenta of the quarks and electrons in 
quark matter that is electrically and color neutral and in 
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where we have simplified p s F upon assuming that M s and 
/i e are small compared to y, by working only to linear 
order in \i e and M^. To this order, electric neutrality 
requires 
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We see from Q that to leading order in M| and fj, e , the 
effect of the strange quark mass on unpaired quark mat- 
ter is as if instead one reduced the strange quark chemical 
potential by M^/(2/i). We shall make this approximation 
throughout. The corrections to this approximation in an 
NJL analysis of a two-flavor crystalline color supercon- 
ductor have been evaluated and found to be small |33l |. 
and we expect the same to be true here. Upon making 
this assumption, we need no longer be careful about the 
distinction between pf's and //'s, as we can simply think 
of the three flavors of quarks as if they have chemical 
potentials 



with 



Md = M« + 2 <5m3 

M« = Pf 

Ms = M« - 25/Z2 



M 2 

<$M3 = <$M2 = = SfJ, 
8/i 



(7) 



(8) 



where the choice of subscripts indicates that 2S[i2 is the 
splitting between the Fermi surfaces for quarks 1 and 3 
and 2(5/i3 is that between the Fermi surfaces for quarks 

1 and 2, identifying u,d,s with 1,2,3. (The prefactor 

2 in the equations defining the Sfx's is chosen to agree 
with the notation used in the analysis of crystalline color 
superconductivity in a two flavor model |29(, in which 
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the two Fermi surfaces were denoted by /i ± Sfi meaning 
that they were separated by 2<5u.) 

As described in Refs. 0,0,11101, BCS pairing intro- 
duces qualitative changes into the analysis of neutrality. 
For example, in the CFL phase fj, e = and fig is nonzero 
and of order M^/fj,. This arises because the construc- 
tion of a phase in which BCS pairing occurs between 
fermions whose Fermi surface would be split in the ab- 
sence of pairing can be described as follows. First, ad- 
just the Fermi surfaces of those fermions that pair to 
make them equal. This costs a free energy price of order 
&\i 2 \i 2 . And, it changes the relation between the chem- 
ical potentials and the particle numbers, meaning that 
the /i's required for neutrality can change qualitatively 
as happens in the CFL example. Second, pair. This 
yields a free energy benefit of order Ag/i 2 , where Ao is 
the gap parameter describing the BCS pairing. Hence, 
BCS pairing will only occur if the attraction between 
the fermions is large enough that A ft 5fi. In the CFL 
context, in which (ud), (us) and (ds) pairing is fight- 
ing against the splitting between the d, u and s Fermi 
surfaces described above, it turns out that CFL pairing 
can occur if A > 4(5/_i = M 2 /(2fi) [ll|, a criterion that 
is reduced somewhat by kaon condensation which acts 
to stabilize CFL pairing |4l"|. In this paper we are con- 
sidering quark matter at densities that are low enough 
(p < M 2 /(2A )) that CFL pairing is not possible. The 
gap parameter Ao that would characterize the CFL phase 
if M 2 and S/j, were zero is nevertheless an important scale 
in our problem, as it quantifies the strength of the attrac- 
tion between quarks. Estimates of the magnitude of Ao 
are typically in the tens of MeV, perhaps as large as 100 
MeV 0. We shall treat A as a parameter, and quote 
results for A = 25 MeV. 

Crystalline color superconductivity can be thought of 
as the answer to the question: "Is there a way to pair 
quarks at differing Fermi surfaces without first equaliz- 
ing their Fermi momenta, given that doing so exacts a 
cost?" The answer is "Yes, but it requires Cooper pairs 
with nonzero total momentum." Ordinary BCS pairing 
pairs quarks with momenta p and — p, meaning that if 
the Fermi surfaces are split at most one member of a 
pair can be at its Fermi surface. In the crystalline color 
superconducting phase, pairs with total momentum 2q 
condense, meaning that one member of the pair has mo- 
mentum p+q and the other has momentum — p + q for 
some p |29l |39j . Suppose for a moment that only u and 
d quarks pair, making the analyses of a two-flavor model 
found in Refs. H OmM IH IM H |M HI (and really 
going back to Ref. j3||) valid. We sketch the results of 
this analysis in the following two paragraphs. 

The simplest "crystalline" phase is one in which only 
pairs with a single q condense, yielding a condensate 

(%p u (x)C~i b il) d (x)) oc Aexp(-2zq • x) (9) 

that is modulated in space like a plane wave. Assum- 
ing that A <c 5^, <C /i, the energetically favored value 
of |q| = q turns out to be q — rjS/j., where the propor- 



tionality constant r\ is given by 77 = 1.1997 [291 139^ . If 77 
were 1 , then the only choice of p for which a Cooper pair 
with momenta (—p + q, p + q) would describe two quarks 
each on their respective Fermi surfaces would correspond 
to a quark on the north pole of one Fermi surface and 
a quark on the south pole of the other. Instead, with 
rj > 1, the quarks on each Fermi surface that can pair 
lie on one ring on each Fermi surface, the rings having 
opening angle ipo — 2cos _1 (l/?7) = 67.1°. The ener- 
getic calculation that determines rj can be thought of as 
balancing the gain in pairing energy as 77 is increased be- 
yond 1, allowing quarks on larger rings to pair, against 
the kinetic energy cost of Cooper pairs with greater to- 
tal momentum. If the A/Sfi — > Ginzburg-Landau limit 
is not assumed, the pairing rings change from circular 
lines on the Fermi surfaces into ribbons of thickness ~ A 
and angular extent ~ A/<5/x. The condensate © carries 
a current, which is balanced by a counter-propagating 
current carried by the unpaired quarks near their Fermi 
surfaces that are not in the pairing ribbons. Hence, the 
state carries no net current. 

After solving a gap equation for A and then evaluating 
the free energy of the phase with condensate © , one finds 
that this simplest "crystalline" phase is favored over two- 
flavor quark matter with either no pairing or BCS pairing 
only within a narrow window 

0.707 < ^ < 0.754 . (10) 
A 

At the upper boundary of this window, A — > and 
one finds a second order phase transition between the 
crystalline and unpaired phases. At the lower bound- 
ary, there is a first order transition between the crys- 
talline and BCS paired phases. The crystalline phase 
persists in the weak coupling limit only if Sfi/Ag is held 
fixed, within the window (|10|) . while the standard weak- 
coupling limit Ao //i — > is taken. Looking ahead to our 
context, and recalling that in three-flavor quark matter 
5fi = Mg/(8/i), we see that at high densities one finds 
the CFL phase (which is the three-flavor quark matter 
BCS phase) and in some window of lower densities one 
finds a crystalline phase. In the vicinity of the second 
order transition, where A — > and in particular where 
A/ 5 fx — > and, consequently given JTDJl, A/A — > a 
Ginzburg-Landau expansion of the free energy order by 
order in powers of A is controlled. Analysis within an 
NJL model shows that the results for A(8/j,) become ac- 
curate in the limit Sfi — > 0.754Ao where A — > 0, as must 
be the case, and show that the Ginzburg-Landau approx- 
imation underestimates A(5/j,) at all S/j, p9ll3r|. 

The Ginzburg-Landau analysis can then be applied 
to more complicated crystal structures in which Cooper 
pairs with several different q's, all with the same length 
but pointing in different directions, arise j3||. This anal- 
ysis indicates that a face-centered cubic structure con- 
structed as the sum of eight plane waves with q's point- 
ing at the corners of a cube is favored, but it does not 
permit a quantitative evaluation of A(8fi) in this case 
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because it predicts a strong first order phase transition 
between the crystalline and unpaired phase, at some Sfi 
significantly larger than 0.754Ao, meaning that there is 
no value of <5/i at which the Ginzburg-Landau approxi- 
mation is under control |36j . 

Our purpose in this paper is to analyze a three-flavor 
analogue of the simplest "crystalline" condensate 
with and without making the Ginzburg-Landau approxi- 
mation A -C Ao. We shall assume weak coupling, namely 
Ao <§; fx, throughout. 

The analysis of neutrality in three-flavor quark mat- 
ter in a crystalline color superconducting phase is very 
simple in the Ginzburg-Landau limit in which A <C 5/i: 
because the construction of this phase does not involve 
rearranging any Fermi momenta prior to pairing, and 
because the assumption A < i5/j implies that the pair- 
ing does not significantly change any number densities, 
neutrality is achieved with the same chemical potentials 
/i e = M 2 / {Aji) and /13 = ^8 = as in unpaired quark 
matter, and with Fermi momenta given in Eqs. Q), (0, 
and 10 as in unpaired quark matter. This result is cor- 
rect only in the Ginzburg-Landau limit, but we expect 
that deviations from it in our NJL analysis are small, 
and so we shall fix the chemical potentials as in unpaired 
quark matter throughout this paper. 

We consider a condensate of the form, 

3 

(iP ai (x)Cj 5 ij 0j (x)) cx J2 A/ e- 2i ^ T e Ia0 e Hj , (11) 

1=1 

where qi, q 2 and q 3 and Ai, A 2 and A 3 are the 
wave vectors and gap parameters describing pairing be- 
tween the (d, s), (u, s) and (u,d) quarks respectively, 
whose Fermi momenta are split by 26 fix, 2b\i2 and 2<5/i3 
respectively. From Q, we see that (5/i 2 = Sfi^ = 
Snx/2 = M^/(8/j). The condensate (jTT|) . first analyzed 
in Ref. [^3, is the natural generalization of the conden- 
sate © to the three-flavor setting, and the natural gen- 
eralization of the CFL condensate (obtained by setting 
qi = q2 = q3 = 0) to the crystalline color superconduc- 
tor setting. 

We shall make the further simplifying assumption that 
Ai = 0. Given that S[ii is twice 6/12 or 8^,3, it seems 
reasonable that Ai <C A2, A3. We leave a quantitative 
investigation of the size of Ai to future work, however. 
With Ai set to zero, the symmetry of the problem is such 
that we expect and find solutions with A2 = A3 = A. 
It should be noted, however, that these simplifications 
are rigorous only in the Ginzburg-Landau limit in which 
A <C S/i. A full investigation of cases in which A ~ <5/i 
requires investigating the Ginzburg-Landau results = 
= 0, [i e = M^/(4/z), Ax = 0, and A 2 = A 3 = A anew, 
as beyond the Ginzburg-Landau approximation these are 
all assumptions, not results. 

In an NJL model in which the quarks interact via a 
point-like four fermion interaction, analyzed at the mean 
field level, if we assume that the only condensate is that 
in Eq. i|ll|) with Ai = the interaction term added to 



the Lagrangian Q can be written simply as 

3 

C A = AJe !2qi ' r e^ e/l ^ M C 75 xj> Pj + h.c. . (12) 
1=2 

Note that since the direction of one of the two wave vec- 
tors is arbitrary, the quantities that have to be deter- 
mined by minimizing the free-energy are A, the magni- 
tude of the two wave vectors, and the angle (f> between 
q 2 and CJ3. We shall see in the next section that, as 
in the two-flavor model, the magnitude of the wave vec- 
tors is given in the Ginzburg-Landau approximation by 
I q 2 1 = I q 3 1 = i]S[i with 77 = 1.1997. And, we shall 
see that, again as in the two- flavor model, A — > at 
Sfj, -> 0.754A , corresponding to M 2 //x = 6.032A . In 
the vicinity of this second order phase transition, the 
Ginzburg-Landau approximation is controlled, as in the 
two-flavor model. However, the angle between the wave 
vectors is a degree of freedom present here that has no 
analogue in the two- flavor model with condensate 10. 
We shall determine it within the Ginzburg-Landau ap- 
proximation in the next section, and check this result in 
Section IV. 



III. GINZBURG-LANDAU ANALYSIS 

The free energy of the crystalline color superconduct- 
ing phase can only depend on the gaps A/ through |A/| 2 , 
because the Lagrangian has U(l) symmetries associated 
with u-, d- and s-quark number conservation. Of course, 
the condensate (|llfl breaks these symmetries, but it does 
so spontaneously. Hence, when the free energy for the 
three-flavor crystalline color superconductor is expanded 
in powers of the A/s, the only terms that can arise up to 
fourth order take the form |A/| 2 , |A 7 | 4 , and |A/| 2 |Aj| 2 
for I 7^ J. As we are setting Ai = 0, this means that 
the Ginzburg-Landau expansion of the free energy can 
be written 

n = a i + |(A 2 +A 2 ) + ^(A 4 +A 4 ) + ^A 2 A 2 +0(A 6 ) , 

(13) 

where we have dropped the absolute value signs as hence- 
forth we will assume that the U(l) symmetries are bro- 
ken such as to give real A's. The contribution Q n is 
the free energy for unpaired ( "normal" ) neutral quark 
matter with ^3 = /is = and [i e = M 2 / (4/x) . Re- 
call that A 2 describes up-strange pairs with momen- 
tum 2q2, A 3 describes up-down pairs with momentum 
2q3, and the three Fermi surfaces are ordered d then u 
then s from biggest to smallest, each separated from the 
next by 28 [i = M 2 /(4/i). The Ginzburg-Landau expan- 
sion is controlled in the limit in which A 2 /i5/i — > and 
A^/Sfi — > 0. In the present context, it was first analyzed 
in Ref. [3l}. The coefficients in the Ginzburg-Landau 
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expansion (|13f) are given by the expressions 
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where q = |q|. Note that in l|13|) the coefficients a and 
/? multiply terms involving only A 2 or only A3. Since 
A2 and A3 separately each describe pairing between two 
fermions with Fermi momenta separated by 2<5/i, this 
means that the expressions for a and (3 are exactly as 
derived in the two- flavor analysis of Ref . [3(| . The coef- 
ficient /3 23 , derived in Ref. j32|, is instead characteristic 
of the 3-flavor case and represents the interaction term 
between the two different condensates. This is the only 
term in which a dependence of the free energy on the 
angle <j> between q 2 and q3 can arise. Note that we have 
written explicitly the it prescription needed to evaluate 
the integral expressions for both (3 and /3 2 3, which we 
have determined by following carefully the derivation of 
(|13|1 from the Ginzburg-Landau expansion of the Green's 
function in the mean field approximation j3(|. The /3 23 - 
integrand in Eq. I|14|) has singular points, making locat- 
ing its poles correctly in the complex plane crucial for 
the evaluation of P23 ■ 

To quadratic order, the free energy (|13|) includes no 
interaction between the condensates A 2 and A3, and so 
is independent of 4>. The expression for a can be analyzed 
as in Ref. [36$ . As long as (3 and /3 2 3 are positive, as we 
shall see below is the case, there is a second order phase 
transition from unpaired quark matter to a phase with 
nonzero A 2 and A 3 at the largest value of <5/i for which 
a = for some value of q. For larger 5 fx, a > for 
all q and the unpaired phase is stable. The transition 
occurs at Sfi = 0.754Ao, where modes with q = for 
r/ = 1.1997 become unstable to condensation, and A 2 
and A3 become nonzero. At lower values of 5fi, modes 
in a band of q have a < making them unstable, but 
the mode with q = r/Sfj, is the one with the most negative 
a and we therefore assume that the condensate involves 
only the modes with q 2 = 173 = r)5[i. 

To evaluate the magnitude of A 2 and A 3 below the sec- 
ond order phase transition where they become nonzero, 
we need the quartic coefficients (3 and /3 2 3. We see that al- 
though the location of the second order phase transition 
is independent of (f>, the magnitude of the condensates 
below the transition can depend on 4> through /3 23 . It is 
convenient to factor out the dependence of (3 and P23 on 
5/i and /x. For q = 1.1997<5/x, then, 



(3 = 2.276 



7T 2 «5/i 2 



(15) 



a 
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FIG. 1: I(cj>), defined in Eq. 1171 . versus the angle 
the wave vectors q2 and q3. 
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between 



and, similarly, we write 



ft: 



23 



2 M 2 

TT 2 6fI 2 



(16) 



defining 



-T(q2,q3) = 5ie 



dp 



-1 



4-7T (it — r?p.q 2 — l)(ie 



?/p.q 3 + 1) ' 

(17) 



From rotational invariance, /(q 2 , 0.3) depends only on 0, 
the angle between q 2 and CJ3. I(4>) can be evaluated nu- 
merically and the result is plotted in Fig^ We note that 
I(4>) is an increasing function of 4> and diverges at 4> = tt, 
i.e. when the two q vectors are antiparallel. The mini- 
mum value of I occurs at <j> — 0, for parallel q vectors. 

We now look for minima of the free energy with A 2 = 
A3 = A, where the free energy can then be written as 



Q = n n + aA 2 + l((3 + p 23 )A 4 



(18) 



The expression for a in Eq. i|14|) yields the result that 
for q = rjSfx, a < for Sfi < 0.7544A , corresponding 
to M 2 /M < 6.035A ~ 151 McV, where we have taken 
Ao = 25 MeV. For values of 5fi where a is negative, the 
solution is found at 



with 



A 



q = n r 



ft + ft 23 



2(/3 + /3 23 ) 



(19) 



(20) 



These expressions can be evaluated using a, [3 and P23 
from Eqs. I] 14(1 . i|15|) and (|16fl respectively. Since (3 + 
/3 2 3 is positive for all choices of the angle <j), we see that 





FIG. 2: Sketch showing where on the Fermi surfaces pairing occurs for condensates in which q2 and q3 are parallel (left 
panel) or antiparallel (right panel). The dark grey (red online) rings on the d and u Fermi surfaces indicate those quarks 
that contribute the most to the (ud) condensate with gap parameter A3 and wave vector q3, which points upward in both 
panels. The light grey (teal online) rings on the u and s Fermi surfaces indicate those quarks that contribute most to the (us) 
condensate with gap parameter A2 and wave vector q2, which points upward in the left panel and downward in the right panel. 
For illustrative purposes, we have greatly exaggerated the splitting between the Fermi surfaces, Sfi/n relative to the values that 
we shall employ in the calculations reported in Section V. 



the transition from the unpaired state to the crystalline 
phase with the simple "crystal" structure that we have 
employed is second order for all choices of <j>. The best 
choice of 4> is the one that minimizes f2, meaning the 
choice with the smallest /?23- From Fig. ^ we see that this 
corresponds to <fi = 0, with q 2 || q3. We shall provide 
plots of A and f2 versus -M 2 //x for various values of 
in Section V, where we shall compare these results to 
results obtained without making the Ginzburg-Landau 
approximation. 

The free energy (|13f) can also be used to analyze the 
free energy of a two-flavor crystalline phase with a single- 
plane wave "crystal structure" in three flavor quark mat- 
ter. Setting A2 = and A3 = A (or equivalently setting 
A3 = and A2 = A) we find a solution with A 2 = \a\/0 
andfi— n„ = — a 2 /(4/3). Like the solution with A2 = A3, 
this solution is neutral in the Ginzburg-Landau limit. 
The solution with A2 = A3 has a lower free energy than 
that with A 2 = if (J3 + P23) < 2/3, meaning #23 < 0. 
One can generalize to other solutions, and rule them 
out, by writing (Aa, A3) as (A r sin 6*, A r cos (9) and then 
rewriting Eq. (|13|) as 



o = o r 



-At 







A 



023-0 



A 2 sin 2 (26>) .(21) 



As long as a < 0, there is a minimum with 9 = ir/A 
(namely A 2 = A3) if P23 < & n d minima at 9 = and 
9 = tt/2 if P23 > 0- From Eqs. iJTHJ and (T^J we sec that 
023 < if I{4>) < 1.138, and from Fig. [I] we see that this 



occurs for < 2.485 radians, or <fi < 142.4°. 

Although the divergence of P23 at = 7r is not phys- 
ically relevant, since for P23 > the lowest free energy 
solution has A2 or A3 vanishing, it is still worth gaining a 
qualitative understanding of the origin of the divergence. 
We see in Fig.|21that there are two pairing rings on the up 
quark Fermi surface, because some up quarks pair with 
down quarks forming Cooper pairs with wave vector 2q3 
and other up quarks pair with strange quarks forming 
Cooper pairs with wave vector 2q2. However, as shown 
in the right panel of Fig. [21 if </> — tt the two pairing 
rings on the up quark Fermi surface are close to coin- 
cident. In the weak-coupling limit in which Sfi/fi — > 
(and Ao — > with Sfi/Ao fixed) these two rings become 
precisely coincident. We attribute the divergence in $23, 
which indicates that antiparallel wave vectors pay an in- 
finite free energy price and hence are forbidden, to the 
coincidence of these two pairing rings. Loosely speak- 
ing, it is as if these up quarks do not know whether to 
pair with their putative strange or down partners and so 
do neither. In contrast, if <f> — as in the left panel of 
the Figure, the two pairing rings on the up Fermi sur- 
face are as far apart as they can be, and (323 and the 
free energy of the state are minimized. This qualitative 
understanding also highlights that it is only in the strict 
Ginzburg-Landau and weak coupling limits that the cost 
of choosing antiparallel wave vectors diverges. If A/S/i is 
small but nonzero, the pairing regions are ribbons on the 
Fermi surfaces instead of lines. And, if Sfi/fi is small but 
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not taken to zero (as of course is the case in Fig. |2J then 
the two ribbons on the up Fermi surface will have slightly 
different diameter, as the Figure indicates. This means 
that we expect that if we do a calculation at small but 
nonzero Ao ~ S[i, and do not make a Ginzburg-Landau 
expansion, we should find some free energetic penalty for 
choosing <fi = it, but not a divergent one. We shall set 
up this calculation in Section IV and see this expectation 
confirmed in Section V. 



IV. NJL ANALYSIS WITHOUT 
GINZBURG-LANDAU APPROXIMATION 

The two-flavor crystalline color superconducting phase 
with a single plane-wave "crystal" structure @ has been 
analyzed in a variety of ways without making a Ginzburg- 
Landau approximation, going back to the work of Fuldc 
and Ferrell (3^. In the QCD context, it was analyzed 
using a variational method in Ref . pflj , using a diagram- 
matic method employing a modification of the Nambu- 
Gorkov formalism in Refs. [13,133, and using the Nambu- 
Gorkov formalism simplified via the High Density Effec- 
tive Theory in Ref. |3lj . 

In the conventional Nambu-Gorkov formalism as ap- 
plied to ordinary BCS pairing, one defines an eight- 
component Nambu-Gorkov spinor 



*(p) 



^(p) 

4> T (- 



(22) 



such that in this basis the pairing between fermions with 
momentum p and — p is described by an off-diagonal 
term in the fermion propagator. The condensate JjJJl, 
however, describes pairing between u quarks with mo- 
mentum p + q and d quarks with momentum — p + q. 
This could be described via a propagator with terms in 
it that are off-diagonal in momentum space, rather than 
merely off-diagonal in "Nambu-Gorkov space" . However, 
it is much easier to change to a basis in which the Nambu- 
Gorkov spinor is written as [30l | 



*(P) = 



/ V>«(p + q) 
Vy(p-q) 
tu(~p - q 

V^!(-p + q) J 



(23) 



The pair condensate is then described by terms in the 
fermion propagator that are off-diagonal in Nambu- 
Gorkov space, occurring in the ipu-4>J and ipd-'tpu en ~ 
tries. In this basis the fermions that pair have p and — p, 
making the propagator diagonal in p-space and the cal- 
culation tractable. One must always keep in mind that it 
is p + q and — p + q that are the momenta of the fermions 
that pair, not p and —p. The variable p is an integra- 
tion variable: in the gap equation or in the expression for 
the free energy, integrating over p sums the contributions 
of all the fermions although of course it turns out that 



only those lying near ribbons on the Fermi surfaces con- 
tribute significantly. Since p is an integration variable, 
we are free to change variables, for example rewriting the 
Nambu-Gorkov spinor as 



*(p) = 



fd{p~ 2q) 
C(-P) 
V^J(-P + 2q) J 



(24) 



The form of the Nambu-Gorkov spinor (|24|l immedi- 
ately suggests that we analyze our three-flavor crystalline 
phase with condensate (|f I \ with Ai set to zero by intro- 
ducing the Nambu-Gorkov spinor 



*(P) 



(Mp) \ 

tpd(p - 2q 3 ) 
fs(p- 2q 2 ) 

C(-P) 
^J(-p + 2q 3 ) 
V^ T (-P + 2q 2 )/ 



(25) 



Furthermore, it also indicates that it will not be possible 
to use this method of calculation if Ai were kept nonzero, 
except for the special case in which qi = q 2 — q 3 . (That 
is, except in this special case which is far from sufficiently 
generic, it will not be possible to choose a Nambu-Gorkov 
basis such that one obtains a propagator that is diago- 
nal in some momentum variable p.) It is thus fortunate 
that, as explained in Section II, it is reasonable on phys- 
ical grounds to begin with Ai = 0, as we do throughout 
this paper. Finally, it seems unlikely that this method 
can be employed to analyze more complicated crystal 
structures analogous to the face-centered cubic structure 
that is favored in the two- flavor case j3||. Investigat- 
ing such possibilities is feasible in the Ginzburg-Landau 
approximation, but the condensate that we are analyz- 
ing, with Ai — and A 2 and A3 each multiplying a 
single plane wave, is the most complex example that 
we currently know how to analyze without making the 
Ginzburg-Landau approximation. 

We now implement the calculation in the basis l|25|) 
using the High Density Effective Theory formalism of 
Refs. [iOIEl' valid in the weak-coupling limit in which 
Ao <C /i. We Fourier decompose the fermionic fields in 
the following nonstandard fashion: 



(26) 

where n is a unit three-vector whose direction is inte- 
grated over, where ki arc three fixed vectors, one for 
each flavor, that we shall specify momentarily and where 
ipia.nix) (resp. ip~ a n( x )) are positive (resp. negative) en- 
ergy projections of the fermionic fields with flavor i and 
color a, as defined in Refs. 0,0. In the usual HDET 
approximation [44|, the vectors are zero and the field 
i/jia,n(x) is used to describe quarks in a patch in momen- 
tum space in the vicinity of momentum p = fin. The 
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introduction of the shift vectors means that now ipia t n(x) 
describes quarks with momenta in a patch in the vicinity 
of momentum /in + k^, with k^ different for different fla- 
vors. To reproduce (|25|l . then, it appears that we should 
choose 



k u = 
kd = 2q 3 
k s = 2q 2 



(27) 



We shall see this choice emerge in a different way below. 

Substituting the expression HU in Eqs. Q and JEJ) 
and neglecting the contribution of antiparticles, the full 
Lagrangian reads 



£ 



dn 

4-7T 



+ (J2 A I e i2t > I - r e^e Ki V n ,iaC75V-n,«e- i ( ki + k J)- r 



1=2 



+ h.c)j 



(28) 



where 6[i%(n) = Pf — fi— k, • n and where the four vectors 
V and V v (the latter used only below) are defined by 
V" = (l,n) and V" = (1,— n). We now see that we 
can get rid of the space dependence in the gap term by 
choosing the shift vectors ki so that they satisfy 



k u 
k, 



k d 

k s 



2q 3 
2q 2 



(29) 



Because the k's were introduced arbitrarily in the decom- 
position i|26|) . the calculation could in principle be done 
with any choice of k's. However, eliminating the space 
dependence in the gap term is an enormous simplifica- 
tion, equivalent to yielding a propagator that is diagonal 
in momentum space, and is what makes the calculation 
tractable. So, we shall always choose k's satisfying (|29ll . 
According to (|29J) . if we choose k M = 0, we recover (12711 . 
However, k u can be chosen arbitrarily as long as k^ and 
k s are then chosen to satisfy (|29|l . This means that the 
choices of k's that get rid of the space dependence in the 
gap term are given by any combination related to 12711 by 
adding any vector to k u and subtracting the same vec- 
tor from kd and k s . The geometric interpretation of two 
examples of choices of k's is described in Fig, [3J As this 
figure illustrates, the freedom to shift k u while keeping 
k u + krf and k u + k s fixed corresponds to the freedom 
to shift integration variable, for example as we did in go- 
ing from (|23|l to (|24|) . In obtaining the results that we 
shall plot in Section V, we shall use the choice l)27|l : how- 
ever, we have checked numerically that different choices 
of k„ with kd and k s satisfying 1|29|) yield the same results 
for the gap parameter and free energy. As Fig. 0] below 
indicates, these different choices yield quite different in- 
termediate stages to the calculation so the fact that we 



find the expected agreement between them is a nontrivial 
check of our numerics. 

We can now employ the Nambu-Gorkov basis defined 
in detail in Ref. [22j given by 



_ J_ ( V>n 



(30) 



where A = 1 ... 9 is a color-flavor index running over the 
nine quarks (three colors; three flavors) and where the 
ip n fields are defined via l|26|l with shift vectors chosen 
as in Eq. I|27|) . In this basis, the full Lagrangian can be 
written in the compact form 



£ 



E 



s a]b( u ) Xb 



(31) 



with 



d a,b — 



(v-£ 



5jJL A (n))5 AB -&AB 

-A AB (V ■£ ~5^ A {~n))5 AB 

(32) 

where £ v = (£o,£\\n) is a four- vector. Here, i\\ is 
the "radial" momentum component of ?/>, parallel to n. 
In HDET, the momentum of a fermion is written as 
(p + f|| )n, with the integration over momentum space 
separated into an angular integration over n and a radial 
integration over —S < £ii < 5. Here, the cutoff S must 
be smaller than /i but must be much larger than Ao, 5fi 
and A. In the results we plot in Section V, we shall take 
li = 500 MeV, 5 = 300 MeV and A = 25 MeV. 

From the Lagrangian H31fl , following a derivation anal- 
ogous to that in Ref. [12], the thermodynamic potential 
per unit volume can be evaluated to be 



O 



~4tt 2 

2A 2 
~G 



E 

£1=1,18 



din 



dn 

-in 



\E*M\ 



12tt 2 



(33) 



where we have set A 2 = A3 = A and where G is a cou- 
pling constant chosen in such a way that Ao = 25 MeV 
in the CFL phase found at /i = 500 MeV with M s = 0. 
In this expression, the E a are the energies of the quasi- 
particles in this phase. They are given by the 18 roots 
of detS" -1 = 0, seen as an equation for £q, with S" -1 the 
Nambu-Gorkov inverse propagator given in l|32(l . The 
quasiparticle energies are functions of £y and n, and also 
depend on the gap parameter A and the wave vectors ki. 
The doubling of degrees of freedom in the Nambu-Gorkov 
formalism means that the 18 roots come in pairs whose 
energies are related by E a (n, £\\) = Ei,(—n,£\\). One set 
of nine roots describes (ipn.djip-n^) and (ipn. Sl tp-n.u) 
pairing, while the other set describes (ipn,u, ip-n,d) and 
("011,11, V'-n.s) pairing. Since n is integrated over, the free 
energy can be evaluated by doing the sum in l|33|) over 
either set of nine roots, instead of over all 18, and multi- 
plying the sum by two. 

In order to determine the lowest free energy state, we 
need to minimize the free energy 17 given in Eq. I|33|l with 
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FIG. 3: Sketches showing how different choices of the shift vectors ki can be interpreted. The left panel shows unshifted d, 
u and s Fermi surfaces, in the absence of pairing. As in Fig. [5] we have exaggerated the magnitude of Sfi/fi for illustrative 
purposes. Now, suppose we have a condensate with q2 = q3 = q, with q a vector pointing upwards. In order to describe the 
(V'n.d, ^-n,u) and (V'n.s, ip-n,u) pairing, we shift the u, d and s Fermi surfaces by — k u , — and — k s respectively (since ip n ,i 
describes quarks in the vicinity of fin + ki) and then reflect (i.e. take n — > — n) the u Fermi surface. The middle panel shows 
the outcome if we follow this procedure with shift vectors given by (12711 . The u Fermi surface is left unshifted (meaning its 
inversion is invisible in the Figure) , and the d and s Fermi surfaces are shifted downwards by 2q. In the right panel, we instead 
choose k u = 2q, which according to 1291 1 then requires k^ = k s = 0. The d and s Fermi surfaces are unshifted. The u Fermi 
surface has been shifted downwards by 2q and then inverted, making it look as if it was shifted upwards. The location of the 
rings on the Fermi surfaces where pairing occurs are determined by the places where the circles cross. In both middle and right 
panels, the pairing rings on the d and s Fermi surface occur at their intersections with the u Fermi surfaces whereas the pairing 
rings on the u Fermi surface are antipodal to where these intersections appear in the Figure, since the u Fermi surface was 
inverted in constructing the Figure. Thus, both the right and middle panels of this figure correspond to the pairing sketched in 
the left panel of Fig. H The same calculation can be done by integrating over the momentum variable of either the middle or 
right panel; the difference between panels is just a change of integration variable. The origins of the momentum variables are 
indicated by the + in each panel. Note that if we had instead chosen to describe the (ip n , u , V'-n.d) an d {ipn,v, ^-n,s) pairing, 
both the middle and right panels would look inverted relative to those given but this difference also corresponds to a change of 
integration variable, in this case n «-» — n. 



respect to the gap parameter A and with respect to </>, the 
angle between q2 and q^. One could also simultaneously 
minimize with respect to fi e , fi^ and fig. And, one could 
allow A2 ^ A3 and minimize with respect to the two gap 
parameters separately. However, in the results that we 
shall present in the next section we shall fix A2 = A3 = 
A, ix e = M^/(4/i) and fiz = fig = 0, as is correct for 
small A and as we have done in the Ginzburg-Landau 
analysis of Section III. 

Before turning to comparing results obtained from the 
calculation presented in this section to those obtained 
with the Ginzburg-Landau approximation developed in 
Section III, we close this section by calculating explicitly 
how the free energy Q of Eq. (|33|) manifests the qual- 
itative features described in Section II and sketched in 
Fig. El with pairing occurring along ribbons of the Fermi 
surfaces. The easiest way to find the regions of momen- 
tum space in which pairing is important is to analyze the 
gap equation, obtained by varying £1 with respect to A. 



This takes the form 

A oc J d£\\ J dn/(n,^||) (34) 

with 

9 n F 

/( n ^ll)=E -g£ Si ^(E a ) ■ (35) 

a—l 

(Either set of nine quasiparticlc energies could be chosen, 
but to make the comparison with Fig. |3] we have used 
those describing (i> n ,d,ip-n,u) an d (t/>a,»,ip-a,u) pairing.) 
In Fig. 21 we plot /, the integrand in the gap equation, 
as a function of t\\ and cos 6* where 9 is the polar angle 
specified by n. (A plot of the integrand in the expression 
for SI in Eq. I|33|) evaluated with A minus that evaluated 
with A = yields a very similar figure.) We have plotted 
/(cos 9, t» ) for two different choices of the shift vectors k; , 
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FIG. 4: Gap equation integrand /, defined in Eq. (I35L as a function of momentum. In each panel, the horizontal axis is 
cosd where 6 is the polar angle specified by n, the vertical axis is In in MeV, and the grey scale indicates the value of /: 
black corresponds to / = and white to the largest values of /. In both panels, /j, = 500 MeV and Sfi — 17.7 MeV, which is 
0.707Ao for Ao = 25 MeV, the interaction strength that we shall use in the next section. In both panels, q2 = q3 = q with 
q — 1.20 Sfi = 21.2 MeV and q pointing in the z-direction. And, in both panels A = 1 MeV. In the left panel, k u = and 
kd = k s = 2q as in the middle panel of Fig. [3] The unshifted u Fermi surface is centered in momentum space, meaning that it 
appears in the left panel as a horizontal line at £m = 0. The shifted d and a Fermi surfaces appear as diagonal lines, with the 
shifted d Fermi surface inside the u Fermi surface at the north pole (cos (9 = 1) and the shifted s Fermi surface outside the u 
Fermi surface at the south pole. In the right panel, k„ = 2q and k<j = k s = as in the right panel of Fig. |3] The unshifted 
d and s Fermi surfaces are centered in momentum space, meaning that they appear in the right panel as horizontal lines at 
l\\ = ±25 ^i. (Note that t\\ is measured relative to where the unshifted u-Fermi surface would have been, shown as a dashed circle 
in the right panel of Fig.|H]that corresponds to in = in the right panel here.) The shifted u Fermi surface appears in the right 
panel here as a diagonal line, outside the d Fermi surface at the north pole and inside the s Fermi surface at the south pole. 
Pairing is most important in the bright white regions, centered where the shifted Fermi surfaces cross. The sketches provided 
in Fig.|5|and particularly in Fig.[3]serve to help visualize the "momentum-space geometry" and pairing regions depicted in the 
present figure. 



corresponding to the middle and right panels in Fig. [3J 
The differences between the two panels of Fig. 0] come 
entirely from the different choices of shift vectors; both 
panels correspond to the same condensate, with q 2 || q3 
as in the left panel of Fig. [2] And, we find excellent 
agreement when we integrate / depicted in either the left 
or the right panel of Fig. 0] to obtain the right-hand side 
of the gap equation (|34[) . and similar agreement when we 
do the integral in Eq. needed to evaluate the free 
energy fl with either choice of shift vectors. 

In both panels of Fig. 21 the bright white pairing re- 
gions near where the shifted and unshifted Fermi surfaces 
cross are clearly visible, as are the jet black blocking re- 
gions near the north and south poles of the Fermi surfaces 
at cos8 = ±1 where no pairing occurs. The pairing re- 
gions are centered at 6 = 67.1°/2 and 6 = 180°-67.1°/2, 
corresponding to cos# = ±0.833. The dark but not black 



regions between Fermi surfaces are places where either 
u-d or u-s pairing is blocked, but the other is allowed. 
Note that even though the formal pairing regions (re- 
gions where / ^ 0) extend far from the Fermi surfaces, 
the bright white regions where the maximal value of / 
is attained are localized near where Fermi surfaces cross. 
So far, all as expected. 

Even with a rather small value of A — in Fig. 0] 
A/5/j, = 0.056 — the ribbons on the Fermi surfaces where 
pairing occurs do not look very narrow. In fact, if we in- 
crease A to the point at which A/dfj, = 0.19, keeping all 
other parameters as in Fig. 21 the blocking regions at the 
north and south poles visible in Fig. ^disappear entirely. 
The pairing "ribbons" become so wide that the ribbon 
centered at 8 = 67.1° expands to encompass 8 — 0, be- 
coming more of a hat for the Fermi surface than a ribbon 
on it. Furthermore, even in Fig. 01 where the pairing rib- 
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M 2 s /\i [MeV] M 2 s /\i [MeV] 



FIG. 5: Plot of A/Ao (left panel) and of the free energy relative to neutral non interacting quark matter (right panel) as a 
function of M 3 / \x for four values of the angle 4> between q2 and q3. The various lines correspond to the calculations done 
in the Ginzburg-Landau approximation described in Section III whereas dots correspond to the NJL calculation of Section 
IV, done without making a Ginzburg-Landau approximation. The full lines (green online) and circles correspond to <j> = 0, 
the dashed-dotted lines (magenta online) and diamonds correspond to (j> = 2tt/3, the dashed lines (red online) and squares 
correspond to (p = 7n/8, the dotted lines (blue online) and triangles correspond to (f> — 31n/32. 



bons on the Fermi surfaces are somewhat narrow, in that 
there are blocking regions at the poles, these blocking 
regions are surrounded by regions of momentum space 
where pairing is quite significant. Just a little distance in 
£|| away from the Fermi surfaces, the angular extent of the 
regions where pairing is significant grows rapidly, becom- 
ing much wider than right at the Fermi surfaces them- 
selves. These are all indications that even though the 
Ginzburg-Landau approximation is formally controlled 
by the parameter A/<5/i, it may break down quantita- 
tively at rather small values of A/Sfi. After all, in the 
Ginzburg-Landau limit A/Sfi — ► the pairing is domi- 
nated by infinitesimally narrow ribbons exactly where the 
shifted Fermi surfaces cross. Using this as a basis for ap- 
proximation cannot yield even a qualitative description 
of the physics once A/8/i ~ 0.2, as even with this small 
a value of A/Sfi the regions where pairing is significant 
are no longer confined to narrow ribbons but have spread 
over considerable regions of the Fermi surfaces. Indeed, 
the extent of the bright white regions where pairing is 
significant in Fig. 0] in which A/5fi = 0.056, indicates 
that the Ginzburg-Landau approximation may cease to 
be quantitatively reliable at values of A/Sfj, below 0.2. 



V. COMPARISONS AND CONCLUSIONS 

In Fig. [S]we compare our results for the gap parameter 
and the free energy in the crystalline color superconduct- 
ing phase calculated in the Ginzburg-Landau approxima- 
tion of Section III with those obtained without making 
this approximation in Section IV. We have done the cal- 
culations with /j = 500 MeV and a coupling strength 



chosen so that Ao, the BCS gap in the CFL phase at 
M s = 0, is 25 MeV. We vary M s , but plot quantities ver- 
sus Ml I M because the most important effect of nonzero 
M s is the splitting between the d, u and s Fermi mo- 
menta given in Eq. J7|), controlled by <5/i = M^/(8/i). 
We analyze the condensate <|11|) with A2 = A3 = A 
and 172 = 93 = q- At each value of M^/fi, we choose 
q = 7/5/1 = T]M^/(&fj,), where 77 = 1.1997. We fix 
M3 = /^s = and /_t e = Mj/(4/j), as appropriate for neu- 
tral three-flavor crystalline quark matter with A/5fi <C 1, 
for which the Ginzburg-Landau approximation is valid. 
We leave investigating the extent to which these chemical 
potentials may shift at larger A to future work. We show 
our results for four values of the angle between q2 and 
q3: = 0, 27r/3, 77r/8 and 3l7r/32. The lines correspond 
to the Ginzburg-Landau analysis of Section III, where we 
have plotted A and Q — 0„ of Eqs. ijTP|) and J50J), using 
Eqs. {13, JTJJ, and JTJJ and using Eq. {T1J to relate a 
to 5fi and hence to M^/fi. The points correspond to the 
NJL calculation of Section IV, where we have minimized 
of Eq. I|33|) with respect to A. 

We see in Fig.[5]that the NJL calculation has a second 
order transition at M^/fi ~ 151 MeV, corresponding to 
5fx ~ 0.754Ao, for all values of the angle cf>. This result is 
in agreement with the Ginzburg-Landau calculation, in 
which the location of the phase transition depends only 
on a, which is independent of <j). We then see that near 
the phase transition, where A/Aq and hence A/Sfj, are 
small, we find good agreement between the NJL calcu- 
lation and the Ginzburg-Landau approximation, as ex- 
pected. For all values of cf>, as A/Aq increases as M^/fi 
is decreased farther from the transition, we see that both 
A and |fl — f2 n | increase more rapidly with decreasing 
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Ml In than predicted by the Ginzburg-Landau calcula- 
tion. When the Ginzburg-Landau approximation breaks 
down, it does so conservatively, undcrprcdicting both A 
and — 0„| for the entire one parameter family of "crys- 
tal structures" parameterized by (f>. (This behavior also 
occurs in the two-flavor model with condensate © 
Furthermore, even where the Ginzburg-Landau approxi- 
mation has broken down quantitatively, it correctly pre- 
dicts the qualitative feature that at all values of / [i the 
most favorable crystal structure is that with = 0. As 
we saw in the discussion of Fig. [2 this can be attributed 
at least qualitatively to the fact that for q2 || q3 the two 
pairing ribbons on the u- Fermi surface are farthest apart. 
We see that the Ginzburg-Landau approximation is use- 
ful as a qualitative guide even where it has broken down 
quantitatively. 

It is evident from Fig.[S]that the extent of the regime in 
which the Ginzburg-Landau approximation is quantita- 
tively reliable is strongly ^-dependent. In the best case, 
which it turns out is <j> = 0, the results of the Ginzburg- 
Landau calculation are in good agreement with those of 
the full NJL calculation as long as A/Aq < 0.25, cor- 
responding to A/dfi < 0.35. On the one hand, this 
looks like a somewhat small value of A/5fi. However, 
it is remarkable that the Ginzburg-Landau approxima- 
tion works so well for this large a value of A/Sfi\ after 
all, we saw in the discussion of Fig. 0] in Section IV that 
for <fi = even with A/Sfi only 0.19 the pairing "ribbons" 
that characterize the Ginzburg-Landau approximation 
have broadened into "pairing hats" encompassing the 
north and south poles of the Fermi surfaces. For larger 
</>, the Ginzburg-Landau approximation yields quantita- 
tively reliable results only for much smaller A. For ex- 
ample, with cj) = 3l7r/32 we have zoomed in on the region 
near the second order phase transition and seen that the 
Ginzburg-Landau calculation does give results in quanti- 
tative agreement with the full NJL calculation, but only 
for A/A < 0.04, corresponding to A/5/J, < 0.05. Why 
does the regime of quantitative validity of the Ginzburg- 
Landau approximation shrink with increasing <ff! One 
possible qualitative explanation is related to the behav- 
ior of the pairing regions displayed for <f> = in Fig. 21 
We saw in that figure that already at quite small A/<5/i, 
pairing is significant over wide swaths of angle in mo- 
mentum space, particularly a little away from the Fermi 
surfaces. Since at </> = tt the two pairing regions on the u 
Fermi surface overlap completely, the closer gets to tt 
the smaller the value of A/5/j, needed for the spreading of 
the pairing regions to become comparable to the separa- 
tion between pairing regions. This qualitative argument 
does suggest that the Ginzburg-Landau approximation 
should break down at smaller A/S/i for larger 0, as we 
find. It does not explain why this breakdown always re- 
sults in A being enhanced relative to that predicted by 
the Ginzburg-Landau approximation. 

A quantitative study of the radius of convergence of 
the Ginzburg-Landau approximation would require eval- 
uating (at least) the A 6 terms, whose coefficients we shall 



generically call 7 . Because we are working in the vicinity 
of a point where a = 0, the first estimator of the radius 
of convergence that we can construct comes by requir- 
ing 7 A 6 < (/3 + P23) A 4 . The results of the comparison in 
Fig.[5]are not conclusive on this point, but they certainly 
indicate that the radius of convergence in A decreases 
with increasing 0, decreasing towards zero for — > tt. 
Given that we see from Fig.^that /?23 (and hence 
increases with and diverges for — > tt, the apparent 
behavior of the radius of convergence in Fig. suggests 
that 7 should diverge even faster for — > tt. 

In Fig. |S|we compare our results for the free energy of 
the three-flavor crystalline color superconductor with the 
most favorable "crystal" structure that we have found, 
namely that with q 2 || q3, to the free energy of the CFL 
and gapless CFL phases [H [H Recall that, as 

discussed in Section I, the gapless CFL phase is known 
to be unstable with respect to counterpropagating cur- 
rents ■ This means that at every value of M% j fj, where 
the CFL phase has been superseded by the gCFL phase 
and where the gCFL phase has lower free energy than un- 
paired quark matter, there must be one or more phases 
with free energy below that of the gCFL phase. We see 
that the crystalline phase we have analyzed, with its par- 
ticularly simple crystal structure consisting of only two 
plane waves, one for the (ud) condensate and one for 
the (us) condensate, has a lower free energy than both 
the gCFL and unpaired phases in a window 128 MeV < 
Mg/fj, < 151 MeV. Recalling that the Ginzburg-Landau 
approximation is reliable over the widest domain pre- 
cisely for the case with q 2 || q3, we note from Fig. [S] that 
the Ginzburg-Landau analysis provides good approxima- 
tion to the results of the full NJL calculation within the 
window 128 MeV < M s 2 /^ < 151 MeV. Hence, our re- 
sult for the extent of this window agrees with that in 
Ref. |37j . which was obtained in the Ginzburg-Landau 
approximation. 

We do not actually expect the condensate we have an- 
alyzed to be the stable, lowest free energy, phase. In- 
stead, we expect that, as in the investigation of two-flavor 
crystalline color superconductivity in Ref. [3(|, crystal 
structures with more plane waves will yield significantly 
lower free energy. It will be interesting to investigate 
which crystal structures are most favorable, and over 
what regime of M^ / (4/i) they have lower free energy than 
the gCFL phase. Looking ahead to such an investigation, 
the implications of our tests of the Ginzburg-Landau ap- 
proximation in the simple case in which we have been 
able to calculate both without and with this approxima- 
tion can be seen as either a glass half empty or a glass 
half full. On the one hand, we find that the approxi- 
mation is quantitatively reliable only for values of A/dfi 
that are small and for some crystal structures (those with 
q 2 and q3 close to antiparallel) very small. This means 
that if there are robust crystal structures with large gaps 
and condensation energies, they may not be described 
reliably within the Ginzburg-Landau approximation. On 
the other hand, we find that even when it breaks down 
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FIG. 6: Plot of the free energy of the three-flavor crystalline color superconducting phase with condensate with Ai = 0, 
A2 = A3, and q2 || q3, compared to the free energy of the CFL and gapless CFL phases. As throughout, we are normalizing 
the interaction to yield a CFL gap parameter at M s = of Ao = 25 MeV. We have plotted the free energy calculated as in 
Section IV, without making a Ginzburg-Landau approximation. We see that this crystalline phase, with its particularly simple 
"crystal" structure, is favored over the other phases shown for 128 MeV < M^/fi < 151 MeV. 



quantitatively the Ginzburg-Landau approximation re- 
mains useful as a qualitative guide, correctly predicting 
that the favored "crystal" structure among our one pa- 
rameter family of possibilities is that with q2 || q3. Fur- 
thermore, the cases in which the domain of reliability of 
the approximation is smallest, those with q2 and q3 close 
to antiparallcl, arc also the cases with the least favorable 
free energy. This suggests that when more complicated 
crystal structures are analyzed, among them the domain 
in which the Ginzburg-Landau approximation is reliable 
may prove wider for those structures which are more fa- 
vorable (i.e. have lower free energy). As an example, 
recall that in our analysis we found narrow domains of 
reliability for the Ginzburg-Landau approximation for <j) 
near it, the same regime where we found that the conden- 
sate we analyzed was anyway not favored because it had 
a larger free energy than one in which only two flavors of 



quarks pair. Finally, in all the cases where we have been 
able to test it, when the Ginzburg-Landau approxima- 
tion breaks down it does so conservatively, underpredict- 
ing the magnitude of A and the favorability of the free 
energy. We conclude that there is strong motivation to 
explore more complicated and realistic crystal structures 
at least initially using this approximation. 
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